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A  thorough  understanding  of  marine  optics  is  a  mission  of  significant  Naval  relevance.  The 
ability  to  use  optical  measurements  (in  situ  and  remote)to  estimate  light  propagation  and 
phytoplankton  growth  rates  is  a  sigi^cant  objective.  This  objective  requires  a  detailed 
understanding  of  the  nature  and  variability  of  source  and  loss  terms,  and  an  ability  to  model  the 
relevant  rates.  Figure  1  summarizes  dominant  state  variables  and  transformations  occurring  in 
marine  optics,  and  the  contributions  by  the  principal  investigator  in  understanding  these 
processes.  Several  significant  problems  previously  have  not  been  resolved  with  respect  to 
measi^g  or  modeling  the  processes  in  Figure  1.  The  long^ange  objectives  of  the  principal 
investigator  is  to  continue  studies  of  the  various  components  and  rate  processes  shown  in  Figure 
1  in  an  effort  to  improve  our  understanding  of  marine  optical  properties  and  photo*physiology 
and  ecology  of  phytoplankton. 


Figure  1.  A  conceptual  model  of  optically  important  processes  in  the  ocean.  The  long-range  objectives 
of  our  research  include  a  detailed  understanding  of  the  processes  and  state  variables  specified  in  this 
figure.  The  state  variables  in  the  boxes  (S1-S6)  correspond  to  dominant  sources  of  optical  variability  in 
the  oceans.  The  principal  investigator  has  dedicated  his  research  effort  to  understanding  the  diverse 
upects  of  this  syuem.  Methods  for  estimating  the  magnitude  and  characteristics  of  SI,  S2  and  ^  have 
included  analysis  of  particle  absorption  using  micropbotometry  (Itnrriaga  et  al.,  1988);  maciophotometry 
(^tchell  and  Kiefer,  1984;  Mitchell  et  al.  1984;  Mitchell  and  Kiefer,  1988a;  1988b;  bfitchell,  1990);  in 
situ  optical  profiling  resulting  in  discovery  of  regional  bioKtptical  relationships  (Mitchell  and  Holm- 
Hansen,  1991a;  Mitchell,  1992);  a  novel  method  to  estimate  SI  fiom  moored  radiometers  (Stramsld  et  al. 
1992);  and  aircraft  (Tania  et  al.,  1990)  or  satellite  observations  of  biomass  (Mitchell  et  al.,  1991).  Optical 
based  models  of  phytoplankton  photosynthesis  for  laboratory  cultures  (Kiefer  and  Mitchell,  1983;  Sosik 
and  Mitchell,  1991)  for  Antarctic  phytoplankton  (Mitchell  and  Holm-Hansen,  1991b)  and  for  Arctic 
phytoplankton  (Cota  et  al.,  1992)  have  been  developed  to  model  process  Kl.  Siegel  et  aL  (1989)  used  an 
in  situ  optical  approach  to  estimate  Kl  and  K2-I-KS  for  populations  in  the  North  Pacific  Gyre.  Carder  et 
(1^1)  bave  used  optical  methods  to  partition  Si  ftom  S2  and  S3  for  remote  sensing.  Sedimentation 
rates  of  Antarctic  and  Arctic  phytoplankton  (K2)  have  been  studied  by  Mitchell  and  Holm-Hansen 
(1991b)  and  Wassman  et  al.,  (1990). 
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OBJECTIVES  FOR  GRANT  N00014.J-1071 


The  wocfc  completed  with  support  of  N000-14-J-1071  was  carried  out  to  improve  our  ability  to 
measure  die  light  absorption  by  the  particles  contained  in  the  boxes  labeled  SI,  S2  and  S3  in 
Hgure  1.  The  principal  goal  of  the  work  in  FY89  was  to  develop  and  validate  a  simple,  practical 
and  accurate  procedure  &  detomination  of  peculate  absorption  coefficients  (Sq)  for  samples 
concentrated  onto  glass  fiber  filters  (Quantitative  Hlter  Tech^que,  QFT).  Includra  in  this  effort 
was  an  analysis  of  the  reproducibility,  limitations  and  statistical  precision  of  the  methods  and  a 
theoretical  treatment  of  empirical  re^ts. 

TASKS  COMPLETED 

During  FY89  a  detailed  study  was  initiated  to  evaluate  the  ability  to  determine  ag  for  dilute 
stupensions  of  pigmented  particles  (e.g.  sea  water)  using  the  Qfr.  A  Pokin-Elmer  Lambda  6 
Hi^  Performance  ^ctrophotometer  with  a  dual  gratins^monochromator  and  an  integrating 
sphere  attachment  was  purchased  to  conduct  the  work.  The  QFT,  integrating  ^here  and  opal 
glass  procedures  were  compared  for  various  species  of  phytoplaidcton  in  culture.  The  integrating 
sphere  attachment  allows  direct  determination  of  ag  for  suspensions.  Di^rent  glass  fiber  filter 
types  and  particle  types  were  analyzed  and  empiridal  algorithms  were  determined  for  correction 
of  the  path  length  aplification  factor  (Beta)  which  is  caused  by  light  scattering  in  the  glass  fiber 
medium. 

Experiments  to  determine  the  pathlength  amplification  factor  for  various  filter  types  was 
completed.  Empirical  algorithms  were  determined  and  optimiz«L  The  absolute  transmittance, 
reflectance  and  absorbtance  for  different  samples  and  densities  was  measured  for  comparison 
with  theoretical  models  of  absorption  by  turbid  samples. 

The  dependence  of  Beta  on  sample  density  was  confirmed;  the  source  of  significant  variance  in 
Beta  is  under  review.  A  previous  optimized  algorithm  for  GF/F  filters  (Mitchell  and  Kiefer, 
1988)  yields  Sp  which  are  reproducible  within  15%  and  which  agree  with  the  integrating  sphere 
results  within  lS%.  New  algorithms  for  several  different  types  of  filters  have  been  developed 
and  have  been  published  (Mitchell,  1990). 

We  have  carried  out  a  detailed  evaluation  of  the  generality,  precision  and  accuracy  of  the 
Quantitative  Filter  Technique  ^FT)  for  determination  of  ^sorption  coefficients  of  aquatic 
particulates  (Mitchell,  1990).  ^e  theoretically  predicted  (Mitchell  and  Kiefer,  1988a)  non> 
linear  nature  of  absorption  ^thlength  amplification  for  particles  retained  on  intensely  scattering 
filters  was  confirmed  and  a  new,  more  direct  algorithm  for  correction  was  developed.  Briefly, 
the  method  was  shown  to  have  (i)  general  utility  such  that  raw  absorbance  determined  on 
different  spectrophotometers  is  comparable;  (ii)  a  precision  of  better  than  ±  10%  for 
determinations  mi  replicates;  and  (iii)  an  accuracy  of  better  than  +  15%  foralgorithms 
determined  for  a  spewed  lot  of  foters.  However,  variability  in  different  filter  lots  of  the  same 
filter  type  can  cause  an  additional  error  of  +  10*15%  if  correction  algorithms  are  not  determined 
for  the  lot  of  filters  used.  ~ 
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Absorption,  fluorescence,  and  quantum  yield  for  growth  in 
nitrogen-limited  Dunaliella  tertiolecta 

Heidi  M.  Sosik  and  B.  Greg  Mitchell 

Marine  Research  Division.  0218.  Scripps  Institution  of  Oceanography.  L  CSD.  La  Jolla.  California  92093 
Abstract 

The  effects  of  steady  state  nitrogen  limitation  on  the  optical  properties  of  Dunaliella  tertiolecta 
were  investigated.  Growih  rate  was  vaned  in  a  continuous  culture  under  constant  irradiance  and 
temperature  with  absorption,  fluorescence,  and  cellular  characteristics  including  pigment  contents 
determined  at  each  steady  state.  The  cellular  concentration  of  Chi  a  increased  with  growih  rate 
while  Chl-a-specific  absorption  and  fluorescence  both  decreased.  In  addition,  the  quantum  yield 
for  growih  varied  by  more  than  a  factor  of  3.5  over  the  growih  conditions  examined,  with  the 
highest  yield  in  the  most  rapidly  growing  cells.  The  decrease  in  magnitude  of  Chl-u-specific  optical 
properties  is  caused  by  pigment  package  effects  and  changes  in  the  abundance  of  accessory  pigments 
relative  to  Chi  a.  Changes  in  absorption  and  fluorescence  properties  are  consistent  with  theoretical 
predictions  for  discrete  packages.  In  addition,  pigment-based  reconstruction  techniques  overesti¬ 
mate  the  magnitude  of  both  in  vivo  and  unpackaged  absorption  by  D.  tertiolecta  under  these  growth 
conditions.  The  observed  variability  in  Chl-a-specific  absorption  and  quantum  yield  is  a  funda¬ 
mental  aspect  of  phytoplankton  physiology  that  should  be  incorporated  into  models  of  oceanic 
primary  production. 


Optical  measurements  from  in  situ  pro¬ 
filing  instruments,  moorings,  aircraft,  and 
satellites  provide  information  at  temporal 
and  spatial  scales  not  possible  with  tradi¬ 
tional  biological  oceanographic  techniques. 
The  potential  for  this  information  to  en¬ 
hance  understanding  of  primary  productiv¬ 
ity  in  the  ocean  depends  on  models  that 
relate  optical  properties  to  photosynthesis. 
Althou^  bio-optical  models  for  primary 
production  and  phytoplankton  growth  have 
been  developed  and  are  being  evaluated 
(Bannister  1979;  Kiefer  and  Mitchell  1983; 
Smith  et  al.  1989),  an  increased  knowledge 
of  the  sources  and  magnitude  of  variability 
in  the  parameters  of  these  models  is  critical 
for  accurate  assessment  of  modeling  capa¬ 
bilities. 

An  important  component  of  these  models 
involves  prediction  of  photosynthetic  pro¬ 
duction  or  growth  from  light  and  chloro¬ 
phyll  data.  Variability  in  this  relationship 
may  be  caused  by  changes  in  Chi  a-specific 
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absorption  or  changes  in  the  quantum  yield 
of  photosynthesis.  The  spectral  Chl-specific 
absorption  coefficient,  ap),*(\)  (units  given 
in  list  of  notation),  is  an  important  param¬ 
eter  of  the  models;  it  provides  an  estimate 
of  the  amount  of  light  absorbed  by  phyto¬ 
plankton  cells  from  knowledge  of  ambient 
irradiance  and  Chi  a  concentration.  The 
photosynthetic  quantum  yield.  is  the 
efficiency  with  which  the  absorbed  light  is 
used  for  C  fixation.  Knowledge  of  the  vari¬ 
ability  in  these  physiological  parameters 
with  growth  conditions  is  important  in 
accurately  applying  bio-optical  models  to 
diverse  oceanographic  regions  and  condi¬ 
tions. 

Modeling  primary  production  over  large 
scales  also  relies  on  estimates  of  pigment 
data  from  remote  sensing.  The  accuracy  of 
current  algorithms  for  retrieving  phyto¬ 
plankton  pigment  concentrations  from  op¬ 
tical  data  is  limited  by  lack  of  knowledge 
concerning  the  variability  .n  pigment-spe¬ 
cific  light  absorption.  This  variability  in¬ 
cludes  both  changes  in  a„,,*(X)  and  changes 
in  the  relative  imp'  iiance  of  phytoplank¬ 
ton.  detritus,  and  dissolved  material  as  ab¬ 
sorbers  tn  the  water  column  (e.g.  Morel  1 988; 
Mitchell  and  Holm-Hansen  1991).  Thus 
physiological  variability  in  pigment-specific 
absorption  can  potentially  be  important  in 
estimating  the  biomass  of  primary  produc- 
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Notation 


Absorption  coefficient  for  phyto- 

plankton  in  vivo,  after  solubiliza¬ 
tion  in  0.5%  Triton  X-IOO,  and 
for  intracellular  matenal.  tn  ' 

Chi  a-specific  absorption  coefficient 

ar,*(A), 

for  phytoplankton  in  vivo,  after 

a„/(X). 

solubilization  in  0.5%  Triton 

X-IOO.  for  a  solution  of  cell  mate¬ 
rial.  and  for  intracellular  material, 
m-  (mg  Chi  a)  ' 

a*(X) 

Specific  absorption  coefficient  for 
pigment  i,  m-  (mg  pigment  /)  ' 

*.(X) 

Backscatter  coefficient,  m  ' 

c 

Intracellular  concentration  of  Chi  a. 
mg  m  ’ 

c, 

Concentration  of  pigment  i  in  a  cell 
suspension,  mg  m  ' 

d 

Cell  diameter,  am 

D(X).  G 

Detection  or  geometric  constant,  di¬ 
mensionless 

£o(X) 

Growth  irradiance.  mol  quanta  m  - 
d  '  nm  ' 

E{\) 

Fluorescence  excitation  irradiance. 
mol  quanta  m  -  s  '  nm  ' 

F 

Fluorescence  from  a  single  cell, 
quanta  s  ' 

F 

Broadband  fluorescence  per  cell,  rel¬ 
ative  quanta  s  '  cell  ' 

F* 

Broadband  fluorescence  per  Chi  a, 
relative  quanta  s  '  (mg  Chi  a)  ' 

a(X) 

Absorption  efficiency,  dimensionless 

a*(X) 

Packaging  parameter,  dimensionless 

X/ 

Wavelength  of  excitation  light  and 
fluoresced  light,  nm 

M 

Specific  growth  rate  (base  e),  d  ' 

Quantum  yield  for  photosynthesis 
and  growth,  mol  C  (mol  quanta)  ' 

Quantum  yield  for  fluorescence,  mol 
quanta  fluoresced  (mol  quanta  ab¬ 
sorbed)'' 

P 

Optical  thickness  along  the  particle 
diameter,  dimensionless 

ers  as  well  as  the  production  realized  by  that 
biomass. 

Chi  (2-specific  absorption  at  the  blue  peak 
and  0ph  have  been  shown  to  vary  in  labo¬ 
ratory  cultures  by  at  least  fourfold  and  ten¬ 
fold,  respectively.  Differences  in  a^*{\) 
among  species  have  been  observed,  and 
variability  within  species  caused  by  changes 
in  growth  irradiance  has  been  well  docu¬ 
mented.  This  variability  in  (2p/,*(X)  is  attrib¬ 
uted  in  part  to  pigment  packaging  effects 
that  vary  with  cell  size  and  pigment  content 
per  cell  (Morel  and  Bricaud  1981;  Mitchell 
and  Kiefer  1 988)  as  well  as  changes  in  the 


abundance  of  accessory  pigments  relative  to 
Chi  a  (Bricaud  et  al.  1983;  Sathyendranath 
et  al.  1987;  Berner  et  al.  1989).  In  addition 
to  variability  in  (2p>,*(X).  there  are  differences 
in  <i>ph  among  species  and  as  a  function  of 
growth  irradiance  for  a  single  species  (Kiefer 
and  Mitchell  1983;  Falkowski  et  al.  1985; 
Sakshaug  et  al.  1989). 

Other  factors  besides  light  can  lead  to 
physiological  adaptations  that  result  in  vari¬ 
ability  in  the  optical  properties  of  phyto- 
planlaon  cells.  Chi  a-specific  absorption  and 
quantum  yield  have  been  shown  to  be  af¬ 
fected  by  N  depletion  in  batch  cultures 
(Welschmeyer  and  Lorenzen  1981;  Cleve¬ 
land  and  Perry  1987).  The  physiological  ef¬ 
fects  of  starvation  may  be  quite  different, 
however,  from  those  observed  under  steady 
state  nutrient  limitation  in  continuous  cul¬ 
ture.  Although  (2p**(X)  can  vary  under  steady 
state  N-limited  growth,  Kiefer  and  Mitchell 
(1983)  hypothesized  that  there  should  be  no 
effect  on  (/>ph.  Recent  studies  of  phytoplank¬ 
ton  in  N-limited  continuous  culture  have 
shown  significant  variability  in  both  api,*(\) 
and  (Aph  (Chalup  and  Laws  1 990;  Herzig  and 
Falkowski  1 989).  Before  the  effects  of  phys¬ 
iological  variability  on  bio-optical  model¬ 
ing  of  primary  production  can  be  evaluated, 
a  more  complete  characterization  of  this 
variability  and  its  sources  is  necessary. 

In  vivo  fluorescence  is  another  optical 
property  widely  used  to  characterize  pri¬ 
mary  producers  in  the  ocean.  Measure¬ 
ments  can  be  made  rapidly  over  a  range  of 
scales.  Accurate  assessment  of  phytoplank¬ 
ton  biomass  is  limited,  however,  by  vari¬ 
ability  in  the  relationship  between  in  vivo 
fluorescence  and  other  variables  like  Chi  a 
concentration.  Previous  studies  of  labora¬ 
tory  cultures  have  documented  that  for  a 
single  species  fluorescence  per  Chi  a  varies 
with  growth  conditions  such  as  light  and 
nutrient  availability  (e.g.  Kiefer  1973; 
Mitchell  and  Kiefer  1 988;  Sosik  et  al.  1 989), 
but  further  investigations  into  the  sources 
of  this  variability  are  necessary  to  reliably 
interpret  fluorescence  signals  from  natural 
samples. 

In  this  work,  we  have  examined  the  effects 
of  steady  state  NO,  -limited  growth  on  pig¬ 
mentation.  absorption,  and  fluorescence  in 
the  marine  chlorophyte  Dunaliella  terti- 
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olecta.  The  roles  of  pigment  packaging  and 
pigment  composition  in  determining  light 
utilization  efficiency  have  been  considered 
and  the  results  placed  in  the  context  of 
theoretical  work  on  absorption  and  fluores¬ 
cence  by  single  cells. 

Materials  and  methods 

Dunaliella  tertiolecta  Butcher  (FCRG 
culture  collection)  was  maintained  in  NOj"- 
limited  continuous  culture  at  four  dilution 
rates.  The  vessel  was  kept  at  22*^;  sterile 
GPM  media  (Loeblich  1 975)  at  full  strength 
except  with  40  /iM  NOs"  was  used  for  di¬ 
lution,  and  the  cultures  were  stirred  contin¬ 
uously  and  bubbled  with  sterile,  water-sat¬ 
urated  air.  Constant  light  was  provided  with 
“cool-white”  fluorescent  lamps  at  an  irra- 
diance  of  165  /imol  quanta  m'-  s  '  for  all 
dilution  rates.  Light  intensity  was  measured 
with  a  Biospherical  Instruments  QSL-100 
quantum  scalar  irradiance  meter.  Absolute 
spectral  irradiance  was  determined  from 
these  measurements  and  the  relative  spec¬ 
trum  for  “cool-white”  lamps  provided  by 
General  Electric.  In  each  case,  analysis  was 
conducted  when  fluorescence  of  the  cell  sus¬ 
pension  remained  constant  for  at  least  2  d 
and  always  after  at  least  10  d  at  a  constant 
dilution  rate.  Specific  growth  rate  m  (base  e) 
was  determined  as  the  dilution  rate  of  the 
culture  (flow  rate  divided  by  the  culture  vol¬ 
ume). 

Subsamples  of  the  cell  suspensions  were 
filtered  onto  Whatman  GF/C  glass-fiber  fil¬ 
ters  and  stored  immediately  at  —  70®C  in 
cryotubes  flushed  with  Nj  gas.  These  sam¬ 
ples  were  later  analyzed  for  pigments  by  re- 
verse-phase,  high-performance  liquid  chro¬ 
matography  (HPLC)  following  the  methods 
described  by  Bidigare  (1989).  Additional 
samples  were  filtered  onto  previously  com¬ 
busted  GF/C  filters  and  frozen  in  combust¬ 
ed  Pyrex  test  tubes.  After  desiccation  of  these 
samples,  carbon  and  nitrogen  contents  were 
measured  with  a  Perkin  Elmer  P-E  2400 
CHN  elemental  analyzer.  Cell  density  and 
equivalent  spherical  cell  volume  were  de¬ 
termined  on  a  Coulter  Electronics  64-chan¬ 
nel,  model  ZH  electronic  particle  counter. 
After  1 5  min  of  dark  incubation,  fluores¬ 
cence  of  the  cell  suspension  was  measured 
on  a  Turner  Designs  model  10  fluorometer. 


Bulk  fluorescence  was  also  measured  with 
a  Sea  Tech  in  situ  fluorometer  after  a  di¬ 
lution  of  ~  1 ;  40  with  filtered  seawater. 

Absorption  properties  were  determined 
with  a  dual-grating  Perkin  Elmer  lambda  6 
UV/VIS  spectrophotometer  equipped  with 
an  integrating  sphere.  Conventional  spec¬ 
troscopic  techniques  do  not  provide  accu¬ 
rate  measurements  of  absorption  for  parti¬ 
cles  in  suspension  (Shibata  1958),  but  the 
integrating  sphere  allows  collection  of  light 
scattered  by  particles  in  the  forward  direc¬ 
tion  thus  improving  absorption  estimates. 
Cell  suspensions  were  concentrated  about 
fivefold  by  centrifugation.  In  vivo  spectral 
absorption,  a^AfA),  was  measured  on  the 
concentrated  samples  in  a  1-cm  cuvette. 
Fresh  culture  medium  was  used  in  the  ref¬ 
erence  beam.  Backscattered  light  is  not  col¬ 
lected  by  the  integrating  sphere,  so  our  es¬ 
timates  are  potentially  greater  than  the  true 
apA(A).  Based  on  the  results  of  Bricaud  et  al. 
(1983),  however,  we  do  not  expect  that  the 
backscatter  coefficient,  b^{K),  was  ever  >  1% 
of  a^A(X).  Absorption  was  also  measured  on 
unconcentrated  cell  suspensions,  but  spec¬ 
tra  were  found  to  be  noisy  for  optical  den¬ 
sities  <0.02.  Comparison  of  the  magnitude 
of  peak  absorption  before  and  after  concen¬ 
tration  show^  a  linear  increase  with  con¬ 
centration  factor,  which  was  calculated  from 
cell  densities  determined  with  the  Coulter 
counter. 

In  addition  to  determining  in  vivo  ab¬ 
sorption,  we  treated  a  subsample  of  the  con¬ 
centrated  susp>ension  with  the  detergent  Tri¬ 
ton  X-100  at  a  final  concentration  of  0.5%, 
with  subsequent  sonication.  This  treatment 
disrupts  the  cells  and  thylakoid  membranes, 
removing  pigment  package  effects  without 
extraction  by  organic  solvents  (Berner  et  al. 
1989).  Spectral  absorption,  ari(X),  was  de¬ 
termined  on  these  solubilized  cells  with  0.5% 
Triton  X-100  in  filtered  seawater  as  a  ref¬ 
erence.  This  treatment  to  remove  package 
effects  should  not  be  considered  a  general 
one  as  we  have  worked  with  other  species 
and  found  variable  success  with  disruption. 
A  further  limitation  of  the  Triton  X  method 
involves  spectral  shifts  in  the  wavelengths 
of  peak  absorption  when  compared  to  in 
vivo  absorption.  We  have  observed  shifts 
toward  shorter  wavelengths  of  1  nm  at  the 
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Table  1.  Cellular  characteristics  of  Dunalielia  tertiolecta  at  several  growth  rates  (base  e)  in  N-limited  con¬ 
tinuous  culture.  Values  in  parentheses  are  I  SD  {n  =  2  for  pigments  and  n  =  3  for  C  and  N  determinations). 
Pigments,  C,  and  N  are  given  in  pg  cell*'. 


Growth  rate  (d 

') 

0.23 

0.61 

1  05 

2.08 

Chi  a 

0.327(0.026) 

0.470(0.017) 

0.730(0.033) 

0.763(0.025) 

Chi* 

0.095(0.010) 

0.153(0.003) 

0.213(0.010) 

0.241(0.004) 

Lutein 

0.149(0.008) 

0.212(0.007) 

0.227(0.017) 

0.281(0.006) 

a-carotene 

0.014(0.002) 

0.0 13(<  0.001) 

0.005(0.002) 

0.017(0.004) 

/3-carotene 

0.096(0.009) 

0.079(0.001) 

0.075(0.018) 

0.105(0.006) 

Neoxanthin 

0.020(0.003) 

0.032(0.001) 

0.042(0.003) 

0.051(0.001) 

Violaxanthin 

0.028(0.001) 

0.037(0.001) 

0.057(0.004) 

0.065(0.001) 

Antheraxanthin 

0.010(0.001) 

0.010(0.001) 

0.012(0.001) 

0.014(0.001) 

C^arbon 

39.68(1.41) 

34.98(0.63) 

26.97(0.42) 

29.68(0.77) 

Nitrogen 

2.07(0.13) 

2.33(0.09) 

2.74(0.16) 

2.72(0.27) 

Cell  diam  (nm) 

6.94(0.06) 

6.91(n.d.) 

6.45(0.15) 

6.96(0.20) 

blue  peak  and  8  am  at  the  red  peak.  To 
compensate  for  these  shifts  when  using 
in  calculations  or  comparison  with 
we  have  shifted  the  spectra  by  1  nm 
for  wavelengths  350-557  nm  and  8  nm  for 
551-750  nm.  The  region  between  551  and 
557  nm  was  repeated  to  prevent  introducing 
a  gap  in  the  spectra;  values  are  relatively 
low  and  uniform  at  these  wavelengths,  so 
this  does  not  distort  the  spectra.  The  spec¬ 
tral  values  of  a7-.r*(^).  however,  are  pre¬ 
sented  without  adjustment  {see  Fig.  2B  and 
Table  2). 

A  subsample  of  the  cell  suspension  was 
also  concentrated  on  Whatman  GF/C  filters 
and  extracted  for  24  h  in  90%  acetone  at 
4*C  in  the  dark.  Spectral  absorption  of  the 
extract  was  measured  as  described  above 
except  that  90%  acetone  was  used  as  the 
reference. 

Results  and  discussion 

Cellular  characteristics— previously 
documented  for  N-limited  growth  of  phy¬ 
toplankton  (e.g.  Caperon  and  Meyer  1972; 
Kolber  et  al.  1988),  we  observed  an  increase 
in  cellular  Chi  a  as  the  growth  rate  of  D. 
tertiolecta  increased  in  the  chemostat  (Table 
1).  This  more  than  twofold  change  in  Chi  a 
was  accompanied  by  a  corresponding  in¬ 
crease  in  Chi  b  per  cell  (Table  1),  resulting 
in  no  trend  in  the  ratio  of  Chi  b  to  Chi  a 
(Fig.  lA).  Increases  were  also  observed  in 
the  absolute  cellular  concentrations  of  the 
carotenoids  (Table  1).  As  found  for  Chi  b. 


the  amounts  of  the  minor  yellow  pigments 
relative  to  Chi  a  changed  very  little  (Fig. 
1 B).  This  constant  stoichiometry  is  not  true, 
however,  for  lutein,  which  decreases  rela¬ 
tive  to  Chi  a  at  higher  u  (Fig.  lA).  Vari¬ 
ability  in  the  ratios  of  a-  and  /9-carotene  to 
Chi  a  was  also  observed  (Fig.  IQ. 

The  increase  in  Chi  a  per  cell  at  higher 
dilution  rate  was  accompanied  by  a  parallel 
increase  in  cellular  N  and  a  decrease  in  cel¬ 
lular  C  (Table  1).  As  previously  observed 
for  this  species  grown  under  N  limitation, 
molar  C :  N  ratios  were  high  and  varied  with 
growth  rate  (Caperon  and  Meyer  1972), 
ranging  from  22  for  p  of  0.23  d"'  to  1 1  at 
IX  of  2.08  d  '.  In  contrast  to  cell  C,  cell  size 
did  not  vary  systematically  with  /x  (Table 
1). 

Absorption  changes— CM  a-specific  ab¬ 
sorption  in  vivo,  ap**(X),  was  highest  at  low 
growth  rales  (Fig.  2A)  that,  as  described 
above,  corresponded  to  decreased  amounts 
of  pigment  per  cell.  This  trend  also  held  for 
absorption  after  treatment  with  Triton  X 
(Fig.  2B)  and  after  extraction  in  acetone  (not 
shown).  Changes  in  ap**(\)  may  be  due  to 
changes  in  the  pigment  package  effect  or  to 
changes  in  the  relative  abundance  of  acces¬ 
sory  pigments  compared  to  Chi  a  (Morel 
and  Bricaud  1981;  Berner  et  al.  1989).  Be¬ 
cause  the  effect  of  packaging  pigments  in  the 
cells  and  chloroplasts  has  b^n  removed  in 
the  Triton  X-treated  samples,  any  changes 
in  flr;r*(X)  should  be  due  only  to  changes  in 
the  pigment  composition  of  the  cells.  At  the 
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Fig.  1.  Abundances  of  accessory  pigments  relative 
to  Chi  0  as  a  function  of  growth  rate  in  N-limited 
Dunaliella  tertioleaa.  Ratios  for  each  pigment  are  shown 
with  symbols  identiAed  in  the  panels.  Standard  devi¬ 
ations  are  smaller  than  the  symbol  size  except  where 
error  bars  are  present. 


Wavelength  (nm) 


Fig.  2.  Chl-o-speciAc  absorption  in  N-limited  Du¬ 
naliella  tertioleaa  at  several  growth  rates.  A.  Absorp¬ 
tion  determined  on  whole-cell  suspensions  showing 
variability  due  to  changes  in  pigmentation  and  pigment 
packaging.  B.  Absorption  of  cells  disrupted  by  treat¬ 
ment  with  Triton  X- 100  so  that  variability  is  due  only 
to  changes  in  the  abundance  of  accessory  pigments 
relative  to  Chi  a. 


blue  absorption  peak,  apA*(X)  is  48%  higher 
at  the  lowest  growth  rate  compared  with  the 
highest.  In  contrast,  the  corresponding  dif¬ 
ference  for  the  Triton  X  absorption  is  only 
24%.  Because  the  relative  variation  in  spe¬ 
cific  absorption  without  package  effects  is 
smaller  than  the  difference  in  vivo,  package 
effects  and  changes  in  pigment  ratios  must 
both  contribute  to  variations  in  ap**(X)  with 
growth  rate.  Our  results  suggest  that,  for  ab¬ 
sorption  of  light  at  the  blue  peak,  about  half 
of  the  observed  variability  in  flp**(X)  is  due 
to  package  effects  with  the  remainder  caused 
by  changes  in  pigment  composition. 
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Changes  in  pigment  composition  are  ex¬ 
pected  to  affect  specific  absorption  at  the 
blue  Chi  a  peak  because  absorption  spectra 
for  accessory  pigments  overlap  the  Chi  a 
spectrum  in  this  region.  Examination  of  the 
pigment  ratio  data  (Fig.  1)  shows  that  re¬ 
sponse  of  the  various  accessory  pigments  to 
increased  N  limitation  differs.  It  is  apparent 
that  the  most  important  changes  leading  to 
increased  aTx*Q^)  at  low  ^  are  in  lutein  and 
/^-carotene  relative  to  Chi  a.  Although  the 
other  pigments  remain  in  fairly  constant 
proportions  to  Chi  a  or  are  present  at  rel¬ 
atively  low  levels,  lutein  and  d-carotene  ra¬ 
tios  are  higher  particularly  at  the  two  lowest 
At.  These  growth  rates  also  correspond  to  the 
highest  arx*Q^)  (Fig-  2B).  In  addition,  at  a 
/i  of  1.05  d"'  the  ratios  are  the  lowest  and 
^tr.Y*(X)  is  also  lowest. 

Although  changes  in  pigment  composi¬ 
tion  play  a  significant  role  in  determining 
in  vivo  absorptive  properties,  as  the  com¬ 
parison  of  absorption  for  whole  cells  and 
solubilized  cells  indicates,  the  packaging  of 
pigments  into  discrete  particles  is  also  im¬ 
portant.  A  theoretical  framework  for  un¬ 
derstanding  this  package  effect  has  been  de¬ 
veloped  (Duysens  1956;  Morel  and  Bricaud 
1981),  and  its  applicability  for  describing 
changes  in  absorption  properties  of  light- 
limited  algal  cultures  has  been  verified 
(Sathyendranath  et  al.  1987;  Bricaud  et  al. 
1988).  Variability  in  ap/.*(X)  observed  for 
nutrient-limited  growth  has  not  previously 
been  interpreted  with  respect  to  the  model 
of  pigment  packaging.  In  addition,  an  im¬ 
portant  parameter  of  the  model,  absorption 
by  intracellular  matter,  acm(X),  has  re¬ 
mained  largely  uncharacterized.  It  is  pos¬ 
sible  to  estimate  the  magnitude  and  vari¬ 
ability  in  flcmfX)  from  Triton  X  absorption 
measurements  and  cellular  pigment  con¬ 
centrations. 

Theoretical  work  on  absorption  by  dis¬ 
crete  particles  allows  absorption  efficiency 
as  a  function  of  wavelength,  Ga(X),  to  be 
predicted  based  on  cell  size  and  a,.„,(X) 
(Morel  and  Bricaud  1981).  The  product  of 
^rY*(X)  and  the  intracellular  concentration 
of  Chi  a  (c)  is  an  estimate  of  ac,„(X), 

a,„(\)  =  arx*W  c.  (1) 

Whereas  Morel  and  Bricaud  (1981)  as¬ 


sumed  a  constant  value  of  a,.„  (=  2  x  lO* 
m“‘)  at  the  blue  peak  for  calculation  of  Q^, 
we  found  that  values  varied  from  1.1  to  2.3 
X  10*  m  ’  mainly  due  to  changes  in  c.  At 
the  red  absorption  peak,  values  were  as  low 
as  0.4  X  10*  m  '  due  to  the  decreased  Chi 
a-specific  absorption  at  red  wavelengths 
compared  to  blue  and  the  lack  of  absorption 
by  accessory  pigments  at  red  wavelengths. 

Morel  and  Bricaud  (1981)  have  described 
absorption  efficiency  as  a  theoretical  func¬ 
tion  of  the  dimensionless  parameter,  p'(X). 
This  parameter  represents  the  optical  thick¬ 
ness  of  a  particle  along  its  diameter  and  is 
defined  as 


p'(X)  =  a,„.(X)  d  (2) 

where  d  is  the  particle  diameter.  The  ab- 
sporption  efficiency,  Qa(\),  is  the  fraction  of 
the  photons  incident  on  the  geometric  cross- 
section  of  the  cell  that  is  absorbed  by  the 
cell.  From  the  absorption  coefficient  for 
phytoplankton  in  a  suspension  ap/,(X),  Qa(X) 
can  be  calculated  as 


<2<,(X)  = 


flpA(X)  F 
(irrfV4)  V 


(3) 


Here  N  is  the  number  of  cells  in  a  volume 
y  of  the  suspension,  and  cells  are  assumed 
to  be  spherical  and  homogenous  (Morel  and 
Bricaud  1981).  The  measurements  that  we 
have  made  allow  independent  calculations 
of  both  p'(X)  and  Qa(X)  for  D.  tertiolecta  at 
each  p  (Fig.  3).  An  index  of  the  package 
effect.  Go*(X),  which  is  the  ratio  of  absorp¬ 
tion  by  discrete  packages  to  the  absorption 
by  the  same  material  dispersed  in  solution, 
has  also  been  presented  by  Morel  and  Bri¬ 
caud  (1981): 

(2/(X)  =  a,**(X)/a,„,-(X).  (4) 

We  applied  this  equation  with  ar.v*(X)  for 
*(X)  in  order  to  calculate  values  of  Q* 
(.v).  For  these  growth  conditions  and  wave¬ 
lengths,  Qa  ranges  from  0.2  to  nearly  0.6, 
while  Qa*  varies  from  almost  0.9  to  as  low 
as  0.5  and  agreement  with  theory  is  excel¬ 
lent  (Fig.  3).  This  variability  suggests  that 
less  severely  N-Iimited  cells  can  absorb  more 
light  and  thus  grow  more  rapidly.  This  in¬ 
crease  in  overall  absorption  by  the  cell  is 
accompanied,  however,  by  a  decrease  in  the 
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Fig.  3.  In  vivo  absorption  efficiency  (Qj  and  the 
packaging  parameter  (Q^*)  at  the  red  and  blue  absorp¬ 
tion  peaks  as  a  function  of  cellular  optical  thickness 
(/)')  for  each  of  the  growth  conditions  examined. 
Q^*,  and  p’  are  estimated  from  measurements  with  Eq. 
1-4  as  described  in  the  text;  the  theoretical  relation¬ 
ships  of  Morel  and  Bricaud  ( 1 98 1 )  are  shown  for  com¬ 
parison  (solid  lines). 


efficiency  with  which  each  molecule  of  pig¬ 
ment  absorbs. 

Pigment  reconstruction  of  absorption— 
Because  of  the  need  for  accurate  estimates 
of  ap/,*(X)  in  bio-optical  models,  a  great  deal 
of  effort  is  currently  being  focused  on  the 
determination  of  this  parameter  for  natural 
samples.  Because  other  absorbing  particles 
are  present  in  these  samples,  direct  mea¬ 
surements  of  particulate  absorption  are  not 
sufficient  to  specify  apA*(X).  This  problem 
has  led  to  the  development  of  alternative 
methods  for  deriving  apA*(X),  including  re¬ 
construction  of  ap/,*(\)  from  HPLC  data  on 
individual  pigment  abundances  (Bidigare  et 
al.  1987).  By  summing  the  contributions  of 
each  of  the  individual  pigments,  ap/.*(X)  can 
be  estimated: 

ap**(X)  =  2  C/Chl  a.  (5) 

Here  a*(\)  is  the  specific  absorption  coef¬ 
ficient  of  pigment  i  at  wavelength  X,  C,  the 
concentration  of  pigment  /  in  the  suspen¬ 
sion,  and  Chi  a  the  concentration  of  Chi  a 
in  the  suspension. 

We  have  made  this  calculation  for  each 
of  the  growth  conditions  studied  using  the 
spectral  specific  absoiption  coefficients  for 
Chi  a,  Chi  b,  and  photoprotective  carot¬ 
enoids  from  Bidigare  et  al.  (1990).  Because 
these  coefficients  are  derived  from  absorp- 
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Fig.  4.  A  comparison  of  several  methods  for  esti¬ 
mating  unpackaged  Chl-a-specific  absorption  for  Du- 
naliella  tertiolecta  grown  at  0.61  d  '.  Absorption  by 
cells  solubilized  in  Triton  X  is  similar  to  that  calculated 
for  a  theoretical  solution  of  cell  material,  but  the  pig¬ 
ment  reconstruction  technique  significantly  overesti¬ 
mates  unpackaged  absoprtion  (details  of  these  methods 
given  in  text).  The  magnitude  of  in  vivo  absorption, 
which  includes  the  manifestation  of  package  effects,  is 
shown  for  comparison.  These  results  are  typical  of  all 
the  growth  rates. 


tion  of  individual  pigments  in  solvents,  any 
in  vivo  package  effects  will  not  be  reflected 
in  flp**(X)  estimated  from  Eq.  5.  For  this 
reason  we  expect  this  estimate  of  flp**(X)  to 
be  comparable  to  the  measured  a7-^*(X).  In 
the  blue-to-green  region  of  the  spectrum, 
reconstructed  ap**(X)  values  are  up  to  2.5 
times  the  true  ap/,*(X)  and  as  much  as  2  times 
aTx*(k),  with  the  largest  differences  at  the 
hipest  growth  rates  (Fig.  4  shows  results  at 
M  =  0.61  d*')-  For  the  various  growth  rates, 
the  mean  value  from  400  to  500  nm  of  the 
reconstructed  apft*(X)  ranges  from  23  to  3 1% 
higher  than  arx*(k). 

As  an  additional  comparison,  we  have 
also  calculated  the  theoretical  absorption  in 
solution  [ajo/*(X)]  following  the  method  of 
Morel  and  Bricaud  (1981).  Inputs  to  this 
calculation  are  the  measured  in  vivo  ab¬ 
sorption  and  the  cell  diameter,  with  the  the¬ 
ory  then  allowing  extrapolation  to  absorp¬ 
tion  at  zero  particle  diameter  (i.e.  for  the 
same  material  in  solution).  For  the  average 
absorption  from  400  to  500  nm,  ajo/*(X)  is 
consistently  lower  than  the  measured 
arx*(k),  but  in  contrast  to  the  reconstructed 
estimate  the  difference  ranges  from  2%  to 
at  most  1 3%. 

The  highest  values  of  ao**(X)  we  estimated 
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with  the  pigment  reconstruction  technique 
correspond  to  values  of  QJS)  that  are  higher 
than  the  theoretical  maximum  of  1.  Un¬ 
doubtedly  package  effects  must  be  included 
before  this  method  can  provide  realistic  es¬ 
timates.  The  reason  for  the  discrepancy  be¬ 
tween  atid  the  estimate  of  ap*Q\) 

from  pigment  reconstruction,  however,  re¬ 
mains  elusive.  The  overestimate  seems  to 
be  due  mainly  to  the  contribution  of  lutein, 
although  the  other  carotenoids  also  play  a 
role.  It  is  possible  that  the  magnitude  of 
absorption  by  these  carotenoids  in  vivo  is 
significantly  different  than  in  oi^nic  sol¬ 
vents  due  to  more  than  package  effects.  Re¬ 
cent  application  of  the  pigment  reconstruc¬ 
tion  technique  in  a  dinofiagellate  includes  a 
correction  for  package  effects  (Nelson  and 
Prezelin  1 990),  but  our  observations  suggest 
that  even  this  modification  would  not  give 
accurate  results  for  D.  tertiolecta.  We  spec¬ 
ulate  that  the  molar  extinction  coefficient 
for  some  of  these  pigments  in  organic  sol¬ 
vents  is  greater  than  the  corresponding  co¬ 
efficient  for  the  same  pigment  in  theoretical 
aqueous  solution.  Whether  this  type  of 
problem  exists  for  other  taxa  with  different 
pigmentation  must  be  tested  with  further 
study  of  laboratory  cultures. 

Fluorescence  changes— Fluorescence  from 
a  single  cell  has  been  described  by  a  theo¬ 
retical  model  that  incorporates  the  role  of 
pigment  packaging  in  modifying  light  ab¬ 
sorption  efficiency  and  intracellular  reab¬ 
sorption  of  fluorescence  (Collins  et  al.  1 985; 
Mitchell  and  Kiefer  1988).  This  model  can 
be  applied  to  conditions  of  broadband  ex¬ 
citation  and  emission  and  can  be  modified 
to  estimate  the  relative  fluorescence  mea¬ 
sured  by  a  specific  detector,  thus  allowing 
the  model  to  be  evaluated  with  our  data  set. 
Collins  et  al.  (1985)  presented  a  functional 
form  for  fluorescence  emission  from  a  cell 
which  is  generalized  here  to  reflect  the  ef¬ 
fects  of  both  excitation  and  emission  wave¬ 
lengths: 

F(\f,  Xy)  =  Xy)  E{\f)  Qa(^e) 

•(irdV4)  <2.*(Xy).  (6) 

Here  0y  is  the  quantum  yield  for  fluores¬ 
cence,  £  is  the  incident  irradiance,  and  X, 
and  Xy  are  the  excitation  and  emission 


wavelengths.  This  generalized  equation 
should  be  applicable  in  cases  where  multiple 
chromaphores  with  different  quantum  yields 
such  as  chlorophyll  and  phycoerythrin  may 
be  contributing  to  the  fluoresced  light  and 
is  valid  for  both  excitation  and  emission 
spectra.  It  is  also  possible  to  express  fluo¬ 
rescence  in  terms  of  ap/,*(X,),  ap/.*(X,),  and 
a,.„,(X,)  if  the  available  information  is  not 
sufficient  to  determine  absorption  efficien¬ 
cies  (Mitchell  and  Kiefer  1988). 

For  spectrally  specified  excitation  irra¬ 
diance  and  absorption  properties,  Ga(X)  and 
Q„*(X),  Eq.  6  can  be  integrated  over  the 
ranges  of  X,.  and  X,  to  give  the  total  flux  of 
fluoresced  quanta  in  a  broad  emission  band 
as  a  result  of  polychromatic  excitation.  In 
order  to  calculate  the  flourescence  from  a 
single  cell  measured  by  a  specific  instru¬ 
ment,  it  is  necessary  to  scale  the  total  flux 
to  include  the  detection  characteristics  D(\,) 
and  the  geometric  constant  G  that  indicates 
the  fraction  of  the  emitted  light  intercepted 
by  the  detector.  The  resulting  expression  is 

F=G  \  I  D(Xy)  F(\„  \r)  dX,  dXy  (7) 

where  F  is  the  fluorescence'  from  a  single 
cell  for  broadband  excitation  and  emission, 
and  F(X„  Xy)  is  as  defined  in  Eq.  6  with  £(X  j 
the  product  of  the  spectral  irradiance  emit¬ 
ted  by  the  excitation  source  and  the  trans¬ 
mission  properties  of  any  filters  or  optical 
components  between  the  source  and  the 
sample.  Similarly,  the  detection  fimction 
Z)(Xy)  is  the  product  of  the  transmission  of 
any  emission  filters  or  optical  components 
and  the  response  of  the  detector  at  each 
wavelength  Xy.  This  fluorescence  can  also  be 
expressed  per  unit  of  Chi  as  F*  which  is 
simply  F  divided  by  the  cellular  concentra¬ 
tion  of  Chi  a. 

Consistent  with  the  foregoing  theory, 
variability  with  growth  rate  documented  for 
specific  absorption  is  also  evident  in  the  flu¬ 
orescence  data.  Despite  differences  in  spec¬ 
tral  characteristics,  results  for  the  Sea  Tech 
fluorometer  (Table  2)  and  trends  observed 
with  the  Turner  Designs  fluorometer  were 
similar.  For  comparison  of  trends,  all  flu¬ 
orescence  data  have  been  expressed  in  rel- 
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Table  2.  Quantum  yield  for  growth  and  optical 
properties  of  N-limited  Dunatiella  certiolecta.  Quan¬ 
tum  yield  has  units  of  mol  C  (mol  quanta)  ';  specific 
absorption  coefficients  are  given  for  the  red  and  blue 
peaks  in  m*  (mg  Chi  a)  '-,  absorption  efficiencies  are 
dimensionless  and  fluorescence  is  expressed  as  relative 
flux  per  cell  (/)  or  per  unit  of  Chi  a  if*). 


Growth  rate  td  ) 

0.23 

0.61 

1.05 

2.08 

(t>ll 

0.013 

0.021 

0.027 

0.048 

a,/(436) 

0.040 

0.037 

0.025 

0.027 

u,/(676) 

0.021 

0.022 

0.018 

0.019 

ar,*(435) 

0.061 

0.056 

0.044 

0.049 

a,f*(668) 

0.024 

0.025 

0.023 

0.025 

(2.(436) 

0.35 

0.47 

0.56 

0.54 

(2.(676) 

0.18 

0.27 

0.40 

0.38 

(2.*(436) 

0.66 

0.67 

0.57 

0.55 

<2.*(676)  0.88 

Sea  Tech  fluorometer 

0.86 

0.78 

0.77 

Measured 

F 

1.00 

1.50 

1.20 

1.40 

F* 

1.00 

1.04 

0.53 

0.60 

Modeled 

F 

1.00 

1.35 

1.32 

1.54 

F* 

1.00 

0.94 

0.59 

0.66 

alive  numbers  scaled  to  be  1  at  the  lowest 
n.  At  the  highest  n,  although  pigment  per 
cell  more  than  doubled  relative  to  the  lowest 
fjL,  Qa  is  higher  by  only  60%  at  the  blue  peak 
and  a  correspondingly  small  increase  is  ob¬ 
served  in  F.  As  pigment  per  cell  increases, 
changes  in  the  pigment  composition  and  in 
the  severity  of  package  effects  are  acting  as 
described  above  to  decrease  ap/,*(X). 

A  similar  decline  in  F*  is  also  evident. 
Because  cell  suspensions  were  kept  optically 
thin  [OE)(675)  <  0.01]  for  all  fluorescence 
measurements,  reabsorption  of  fluoresced 
light  within  the  suspension  cannot  account 
for  the  observed  variability  in  fluorescence. 
For  the  Turner  Designs  measurements,  red 
absorption  by  the  cells  in  suspension  was 
about  equal  to  absorption  by  the  water  me¬ 
dium  and  an  order  of  magnitude  lower  than 
absorption  by  the  water  for  measurements 
with  the  Sea  Tech  fluorometer.  Although  it 
is  possible  that  variability  in  the  quantum 
yield  of  fluorescence  (number  of  photons 
fluoresced  per  number  absorbed)  could  oc¬ 
cur  for  cells  from  different  growth  condi¬ 
tions,  the  trends  we  observed  suggest  that 
variability  in  the  efficiency  with  which  light 
is  absorbed  is  an  important  factor. 


From  observed  values  of  absorption  and 
cell  size  and  the  known  characteristics  of 
the  Sea  Tech  fluorometer,  we  have  also  pre¬ 
dicted  the  relative  fluorescence  with  Eq.  6 
and  7.  Because  we  do  not  have  information 
concerning  variability  in  0,,  as  a  first  ap¬ 
proximation  we  have  assumed  that  it  is  con¬ 
stant  for  the  excitation  and  emission  con¬ 
ditions  used  and  that  it  does  not  vary  over 
the  growth  conditions  examined.  Because 
we  do  not  know  the  absolute  values  of  (t>, 
and  G.  fluorescence  values  have  been  pre¬ 
sented  as  relative  numbers  scaled  to  maxi¬ 
mum  of  1  at  M  of  0.23  d  ’.  As  indicated 
above,  these  calculations  involve  not  only 
spdctral  values  of  the  absorption  factors  for 
the  cells  but  also  the  spectral  output  of  the 
lamp  and  the  spectral  transmission  of  the 
blue  excitation  filter,  the  acrylic  window, 
the  red  emission  filter,  and  the  gold  mirror, 
as  well  as  the  spectral  sensitivity  of  the  de¬ 
tector  (see  Bartz  et  al.  1988). 

Trends  in  values  of  F*  predicted  with  Eq. 
6  and  7  agree  well  with  actual  measurements 
(Table  2,  H  =  0.96),  which  supports  our 
hypothesis  that  changes  in  absorption  effi¬ 
ciency  are  an  important  source  of  variability 
in  fluorescence  per  Chi  a  and  demonstrates 
The  validity  of  the  model  for  broadband  ex¬ 
citation  and  emission  of  fluorescence.  Re¬ 
gression  with  the  measured  and  modeled 
values  as  independent  and  dependent  vari¬ 
ables.  respectively,  reveals  an  intercept  of 
0.20.  With  so  few  points  it  is  not  reasonable 
to  attach  much  significance  to  this  intercept, 
but  it  is  possible  that  it  reflects  a  change  in 
<i)f  with  growth  conditions.  Only  a  larger 
range  of  experimental  conditions  and  more 
detailed  work  on  variability  in  F*  can  test 
this  hypothesis.  Further,  the  broadband  ap¬ 
proach  presented  here  is  somewhat  crude, 
and  spectral  fluorescence  methods  such  as 
those  of  Mitchell  and  Kiefer  (1988)  will  be 
necessary  to  resolve  the  details  in  model 
verification. 

The  variability  in  fluorescence  properties 
observed  must  be  considered  in  the  context 
of  the  increasing  use  of  in  situ  fluorometers 
to  estimate  phytoplankton  biomass.  Al¬ 
though  with  proper  calibration  this  method 
may  provide  useful  information,  no  simple 
relationship  can  be  expected  to  exist  be¬ 
tween  fluorescence  and  biomass  estimates 
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such  as  Chi  a  concentration.  Both  variabil¬ 
ity  with  growth  conditions  for  a  sin^e  spe¬ 
cies  and  the  potential  for  interspecific  dif¬ 
ferences  must  be  considered.  In  addition, 
the  lack  of  knowledge  concerning  variability 
in  the  quantum  yield  for  fluorescence  con¬ 
tinues  to  limit  applications  such  as  the  use 
of  natural  fluorescence  to  estimate  produc¬ 
tivity. 

Quantum  yield  for  growth —The  quantum 
yield  at  the  irradiance  under  which  cells  have 
been  grown,  can  be  determined  from  a 
spectral  version  of  the  model  of  Kiefer  and 
Mitchell  (1983); 

M  =  Chi :  C  J  £„(X)  dX.  (g) 

Here  u  is  the  specific  growth  rate  of  the  cells, 
£o(X)  is  the  spectral  growth  irradiance,  and 
Chi ;  C  is  the  cellular  Chi  a  :  C  ratio.  Solving 
this  equation  for  at  each  of  the  growth 
rates  examined  results  in  nearly  fourfold 
variability  in  with  values  decreasing  as 
N  limitation  increases  (Table  2). 

A  potential  source  of  variability  in  <t>^  or 
<^ph  is  changes  in  ap^*(\)  due  to  pigment- 
ratio  shifts.  If  the  abundance  of  photopro- 
tective  pigments  increases  relative  to  photo¬ 
synthetic  pigments,  apt,*{X)  is  expected  to 
increase  but  the  energy  available  for  C  fix¬ 
ation  would  decrease.  The  result  is  a  de¬ 
crease  in  with  higher  apf,*{X),  consistent 
with  our  observations  of  increased  abun¬ 
dance  of  lutein  and  /3-carotene  relative  to 
Chi  a  at  low  p.  This  increase  in  ap**(X)  can 
only  account  for  ~25%  of  the  change  we 
oberved  in  </>„,  however,  because  varied 
by  more  than  threefold  and  pigment-ratio- 
induced  changes  in  ap**(X)  were  much 
smaller  (see  Table  2). 

Many  factors  besides  ap**(X)  may  con¬ 
tribute  to  a  decrease  in  including  changes 

in  respiration  relative  to  photosynthesis  with 
growth  rate.  Although  we  have  not  mea¬ 
sured  respiration  in  this  study,  it  is  unlikely 
that  increased  respiration  relative  to  growth 
could  account  for  all  of  the  decrease  in 
observed  at  low  u-  For  N-limited  growth  of 
Pavlova  lutheri,  Chalup  and  Laws  (1990) 
found  a  more  than  twofold  increase  in  0ph 
with  p  and  also  demonstrated  an  increase 
in  (/>ph  with  decreased  irradiance  at  the  same 
dilution  rate.  For  N  limitation  of  Isochrysis 


galbana  at  a  single  irradiance,  Herzig  and 
Falkowski  (1989)  observed  a  range  in 
similar  to  our  observations  and  additionally 
found  that  variability  in  0ph  was  85%  of  that 
for 

Additional  sources  of  variability  in  quan¬ 
tum  yield  were  not  investigated  in  this  study. 
Changes  in  energy  transfer  efficiency,  in  the 
reaction  center  turnover  time,  or  in  the 
abundance  of  pigment  in  photosystem  1  rel¬ 
ative  to  photosystem  2  could  all  contribute 
to  changes  in  Although  Cleveland  and 
Perry  (1987)  have  implicated  these  mech¬ 
anisms  to  explain  quantum  yield  changes  in 
transition  from  nutrient-replete  to  nutrient- 
starved  conditions,  results  may  differ  for 
cultures  allowed  to  adapt  to  a  steady  state 
N-limited  growth  condition.  Using  pump 
and  probe  fluorescence  techniques  and  im- 
munochemistry,  Kolber  et  al.  (1988)  sug¬ 
gested  that  N  limitation  in  several  species, 
including  D.  tertiolecta,  leads  to  a  decrease 
in  the  energy  conversion  efficiency  of  photo¬ 
system  2  caused  by  a  loss  of  functional  re¬ 
action  centers  or  inefficient  energy  transfer 
to  the  reaction  centers. 

All  of  these  results  suggest  that  studies 
considering  photosynthetic  efficiency  as  in¬ 
dependent  of  nutrient  availability  must  be 
viewed  with  caution.  Kiefer  and  Mitchell 
(1983)  hypothesized  that  <f>ph  is  a  function 
only  of  irradiance,  consistent  with  their 
model  for  growth  of  Thalassiosira  fluviatilis 
in  which  spectrally  averaged  Op,,*  was  as¬ 
sumed  to  be  constant.  Sakshaug  et  al.  ( 1 989) 
used  a  similar  fitting  technique  to  model 
growth  rate  under  a  large  number  of  light- 
and  nutrient-limited  conditions  and  found 
a  constant  value  of  the  product  ap>,*  -0ph  to 
be  adequate  for  a  given  irradiance.  The  con¬ 
clusion  that  (f>ph  is  only  a  function  of  growth 
irradiance  alone  and  not  of  nutrient  limi¬ 
tation.  however,  has  not  been  substantiated 
by  studies  in  which  light  absorption  and 
photosynthetic  efficiency  have  l»en  mea¬ 
sured  (e.g.  this  study;  Kolber  et  al.  1988; 
Herzig  and  Falkowsfe  1989). 

Conclusion 

It  is  evident  that  variability  in  the  phys¬ 
iological  responses  of  phytoplankton  to 
growth  conditions  cannot  be  ignored  in  bio- 
optical  characterization  of  ocean  waters  both 
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with  respect  to  productivity  and  optical 
properties.  Models  for  the  remote  estima¬ 
tion  of  pigments  or  for  propagation  of  light 
through  the  water  column  rely  on  accurate 
knowledge  of  the  relationships  between  pig¬ 
ments  and  in  situ  optical  properties.  For 
example,  recent  work  suggests  that  predic¬ 
tions  with  current  algorithms  for  pigment 
estimation  may  lead  to  significant  errors  for 
the  Antarctic  Ocean  (Mitchell  and  Holm- 
Hansen  1991).  These  errors  result  from 
variability  in  ap/,*(X)  and  in  the  relative 
abundance  of  phytoplankton  and  other  ab¬ 
sorbing  material  in  the  water  column.  Es¬ 
timating  photosynthetic  biomass  from  in 
situ  fluorescence  measurements  is  equally 
complex  and  can  be  significantly  affected  by 
the  physiological  state  of  the  phytoplank¬ 
ton. 

Bio-optical  models  of  phytoplankton 
growth  and  productivity  must  also  incor¬ 
porate  physiological  variability  if  light  ab¬ 
sorption  is  to  be  accurately  related  to  C  fix¬ 
ation.  Working  with  a  model  similar  to  that 
of  Kiefer  and  Mitchell  (1983),  Sakshaug  et 
al.  (1989)  found  that,  compared  to  other 
model  parameters,  changes  in  the  value  of 
Up/,*  caused  relatively  large  errors  in  mod¬ 
eled  growth  rate.  Increasing  evidence  of  large 
variability  in  the  optical  properties  of  phy¬ 
toplankton  with  growth  conditions  suggests 
that  modification  of  existing  models  to  in¬ 
corporate  physiology  and  environmental 
conditions  besides  light  intensity  must  be  a 
high  priority.  A  major  obstacle  to  the 
achievement  of  this  goal  continues  to  be  a 
lack  of  comprehensive  physiological  studies 
that  include  simultaneous  measurement  of 
pigments  and  optical  properties  for  a  wide 
range  of  species  and  growth  conditions. 
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ABSTRACT 

Particulates  in  aquatic  environments  dominate  variability  in  inherent  and  apparent  optical 
properties.  Dilute  concentrations  have  made  the  quantitative  determination  of  particulate  absorption  a 
relatively  difficult  problem.  A  quantitative  filter  technique  (QFT)  has  been  described  previously 
which  allows  determination  of  the  particulate  absorption  coefficient  (a.(A))  for  samples  concentrated  on 
filters  by  correcting  for  the  pathlength  amplification  effect  (fi).  Details  of  the  generality,  accuracy  and 
precision  of  the  procedure  have  not  been  reported  previously. 

A  spectrophotometer  equipped  with  an  integrating  sphere  accessory  was  used  to  study  the 
optical  density  of  phytoplankton  suspensions  and  colored  polystyrene  beads,  ODj(A),  and  the  optical 
density  of  the  same  suspensions  on  a  variety  of  common  filter  types,  ODf<A).  The  ratio  of 
ODf<A)/ODj(A)  confirms  that  multiple  scattering  in  fitters  leads  to  variable  0,  a  non-linear  function  of 
ODf<A).  Studies  of  ODf(X)  when  the  filters  are  mounted  in  different  positions  in  a  single 
spectrophotometer,  and  for  standard  and  integrating  sphere  spectophotometers  indicate  that  the 
procedures  presented  are  of  general  applicability.  With  proper  care  in  baseline  correction  and  sample 
preparation,  spectra  of  ap(A)  can  be  calculated  from  measurements  of  ODf(A)  for  aquatic  particles  of 
diverse  size  and  refractive  index  with  an  accuracy  of  better  than  ±  15%.  A  comparison  of  filter  types 
commonly  used  for  aquatic  research  indicates  that  different  algorithms  are  required  for  estimation  of 
ap(A)  from  ODf(A)  measurements  using  different  filter  types. 

I.  INTRODUCTION 

The  spectral  absorption  coefficient  (a(A))  is  an  inherent  optical  property  of  aquatic  systems^<^. 
As  an  inherent  optical  property,  a(A)  can  be  partitioned  into  its  subcomponents: 

a(A)  -  ♦  a,(A)  ♦  ap(A)  (I) 

where  a^,  a,  and  ap  are  the  portion  of  the  coefficient  attributable  to  water,  soluble  constituents  and 
particulates,  respectively. 

Variations  in  the  apparent  optical  properties  of  water,  including  the  spectral  diffuse  attenuation 
coefficient  k(A)  and  irradiance  reflectance  R(A),  are  determined  by  variations  in  the  spectral  backscatter 
coefficient  and  afA)^'^: 
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k(A)  D  a(A)  ♦  b^(X) 

R(A)  -  0.33  h^(X)  /  a(A)  (3) 

where  D  is  the  radiance  distribution  function.  Thus,  understanding  variability  in  a(A)  is  crucial  for 
understanding  variability  in  these  (and  other)  apparent  optical  properties  of  aquatic  systems.  Since 
ap(A)  dominates  the  variability  of  a(A)  in  oceanic  waters  with  little  terrestrial  influence,  it  is  among  the 
most  important  parameters  which  must  be  measured  for  a  complete  description  of  the  optical  properties 
of  aquatic  systems. 

The  method  of  concentrating  marine  particles  on  filters  for  measurement  of  optical  density 
(ODf<A))  was  originally  applied  in  field  studies  nearly  thirty  years  ago  by  Charles  Yentsch^.  ODf(A)  is 
calculated  from  the  ratio  of  the  transmittance  through  a  reference,  T^(A),  and  sample  filter,  Tj(A); 

ODf(A)  -  -Logjo  (T5(A)/T^(A))  (4) 

Kiefer  and  Soo  Hoo^  demonstrated  that  the  observed  magnitude  of  ODf(A)  was  greater  than  one  would 
predict  based  on  the  geometric  pathlength  of  the  ratio  of  volume  filtered  (V^)  to  clearance  area  of  the 
filter  (A^).  They  concluded  that  the  optical  pathlength  was  greater  than  the  geometric  pathlength  due 
to  scattering  by  the  filter.  Butler^  defined  the  phenomenon  of  pathlength  amplification  for  absorption 
of  material  in  scattering  suspensions: 


^  •  /q  /  /g  (5) 

where  and  /g  are  the  optical  and  geometric  pathlengths,  respectively.  The  value  of  L  for  particles 
concentrated  on  filters  is  simply  the  ratio  of  the  volume  filtered  to  the  clearance  area  of  the  filter: 

/g  -  Vf  /  (6) 

where  A^  -  tt  (7) 

and  where  r  is  the  radius  of  the  portion  of  the  filter  which  contains  the  particles  removed  from  the 
suspension.  From  the  above  relationships,  the  absorption  coefficient  ap(A)  can  be  calculated: 

ap(A)  -  2.3  OD^A)  /  (/g  P)  (8) 

The  factor  2.3  converts  LogjQ  to  Log^. 

Mitchell  and  Kiefer demonstrated  the  non-linear  relationship  between  0  and  ODf<A),  and 
suggested  that  simple  empirical  relationships  could  be  used  to  make  quantitative  estimates  of  ap(A)  from 
measurements  of  ODf(A).  The  concept  of  this  quantitative  filter  technique  (QFT)  for  correction  of  0  to 
estimate  ap(A)  for  dilute  suspensions  has  been  widely  applied* but  there  has  been  little 
uniformity  with  respect  to  correction  methodology.  A  review  of  inconsistencies  in  methods  presented 
in  the  literature  indicates  that  the  various  methods  might  differ  by  as  much  as  two-fold  in  the  final 
estimate  of  ap(A).  For  example,  Kiefer  and  Soo  Hoo^  based  their  estimates  of  ^  on  a  comparison 
between  ODf(A)  determined  on  field  samples  and  diffuse  absorption  coefficients  determined  for 
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cultures^^.  Subsequently,  others  have  used  the  same  ft  for  estimates  of  the  absorption  coefficient  which 
may  be  as  small  as  half  the  value  of  the  diffuse  absorption  coefficient^.  An  evaluation  of  the  QFT 
methodology,  its  generality,  precision  and  accuracy  is  presented  in  this  report. 

2.  METHODS 

Cultures  of  unicellular  algae  were  grown  in  nutrient  enriched  media  and  the  optical  density  of 
suspensions  (OD,(A))  was  measured  in  a  spectrophotometer.  A  variety  of  cell  sizes  and  taxonomic  types 
were  used  including  chlorophytes  (Dunaliella  tertiolecta^  Chlorella  sp.\  Nemochloroxis  sp.);  a 
bacillariophyte  (Thalasiossira  fluviatilis);  and  a  cyanobacterium  (Synechoeoeia  sp.).  Cell  sizes  ranged 
from  2-20  pm.  Blue  polystyrene  beads  (2.85  pm.  Polysciences,  Inc.)  were  also  used  as  a  particle  type 
with  high  refractive  index  relative  to  water  (n  1.2),  compared  to  the  algae  (n  -  1. OS- 1.1).  The 
refractive  index  of  polystyrene  is  approximately  1.4S  relative  to  air.  In  addition,  a  sample  of  natural 
sea  water  collected  from  the  Scripps  Institution  of  Oceanography  pier  was  concentrated  by  filtration 
and  resuspended  for  measurement  of  ODj(A)  according  to  the  methods  of  ICirk^^. 

Measurements  of  OD(A)  were  made  with  a  dual  grating  double  beam  Perkin-Elmer  Lambda  6 
spectrophotometer  equipped  with  an  integrating  sphere.  Data  acquisition  was  accomplished  via  a  digital 
IEEE  488  interface  between  the  spectrophotometer  and  a  personal  computer.  The  spectrophotometer 
automatically  corrected  for  the  baseline  which  was  stored  prior  to  analyses  of  different  types.  For 
suspensions,  the  baseline  was  a  filtrate  of  the  growth  medium.  For  the  beads,  distilled  water  was  used. 
For  the  filters,  the  baseline  was  determined  for  filters  which  were  saturated  with  the  filtrate  of  the 
growth  medium  or  distilled  water,  in  the  case  of  analysis  of  beads.  Values  of  OD,(A)  were  directly 
determined  for  suspensions  of  cells  or  beads  using  the  integrating  sphere  attachment  with  the 
appropriate  medium  in  the  reference  beam.  For  phytoplankton  and  natural  particulates  from  sea  water, 
spectra  were  set  to  0.0  at  750  nm.  The  suspensions  were  concentrated  onto  various  filter  types  for 
determination  of  ODf(A),  with  the  appropriate  filter  type  in  the  reference  beam.  The  reference  and 
sample  filter  were  kept  saturated  during  the  analyses.  This  was  ensured  by  cutting  a  section  of  1  cm 
width  from  the  center  of  the  filter,  and  placing  the  section  into  a  1  cm  cuvette  with  a  drop  of  the 
appropriate  fluid.  Caution  must  be  exercised  to  maintain  saturation,  but  not  oversaturation  which  can 
result  in  resuspension  of  particles  off  the  filter.  Filter  types  studied  included  Whatman  glass  fiber 
types  GF/F,  GF/C,  934AH  and  Miliepore  cellulose  acetate  type  HA. 

The  relationships  between  ODj(A)  and  ODf(A)  were  studied  using  statistical  curve  fitting 
procedures  in  SYSTAT,  a  commercial  software  package  available  for  personal  computers.  Various 
procedures  to  fit  the  observed  non-linearity  in  the  relationship  between  OD,(A)  and  ODf(A)  were 
investigated,  and  a  simple  quadratic  equation  (second  order  polynomial)  was  found  to  be  the  most 
satisfactory.  Thus,  the  magnitude  of  OD,(A)  may  be  estimated  directly  from  ODf<A)  by  a  relationship 
of  the  form: 


OD,(A)  -  a  ODj<A)  ♦  b  IODj<A)j2 


(9) 
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2.  RESULTS  AND  DISCUSSION 


2.1  Generality  of  ODf(A)  determinations 

Of  paramount  concern  with  respect  to  the  utilization  of  the  QFT  is  the  issue  of  the  generality  of 
determination  of  ODf<A)  using  different  spectrophotometers  or  different  optical  geometry  for 
determining  T^(X)  and  Tj(X).  Provided  that  the  post-sample  path  length  for  the  reference  and  sample 
positions  is  equivalent,  the  primary  difference  between  systems  or  configurations  is  the  angle  of 
acceptance  of  the  scattered  light  transmitted  from  the  diffusing  sample  to  the  detector.  Changing  the 
angle  of  acceptance  from  the  sample  to  the  detector  could  hypothetically  cause  variations  in  the  OD^a) 
observed  for  the  same  sample. 

An  experiment  was  done  to  evaluate  the  effects  of  varying  acceptance  angle  of  the  post-sample 
beam.  ODf<A)  was  determined  for  3  separate  positions  along  the  optical  axis  in  the  sample 
compartment  of  the  Lambda  6  spectrophotometer  configured  without  the  integrating  sphere.  A  fourth 
geometry,  attained  by  determining  ODf(A)  using  the  integrating  sphere,  resulted  in  an  angle  of 
acceptance  of  nearly  180*.  For  each  geometry,  a  baseline  was  determined  by  placing  blank  filters  in 
the  sample  and  reference  positions.  ODf(A)  was  then  determined  for  replicate  samples  which  had  the 
same  volume  of  a  culture  of  D.  lertiolecta  filtered  onto  GF/F  filters.  Results  of  this  experiment  are 
presented  in  Figure  I.  There  was  no  significant  difference  between  values  of  ODf(A)  for  the  different 
positions,  or  for  the  integrating  sphere  determination.  Differences  between  replicate  samples  was 
comparable  to  differences  for  the  various  measurement  geometries. 

2.2  Detailed  studies  of  GF/F  and  GF/C  fitter  types 

The  most  commonly  used  filters  for  either  oceanographic  field  work  or  laboratory  work  are  of 
the  Whatman  GF/F  or  GF/C  type.  Among  the  concerns  raised  by  other  investigators  regarding  the 


Figure  I .  A  comparison  of  observations  of 
ODf(A)  for  different  optical  geometries. 
After  storing  a  baseline  for  each  sample 
position,  ODf(A)  was  determined  for 
duplicate  filters  containing  the  same  volume 
filtered  of  D.  lertiolecta.  Samples  were 
placed  along  the  optical  axis  at  the  front, 
center  and  back  of  the  sample  compartment 
(IS  cm  range)  and  in  the  integrating  sphere 
attachment.  Each  curve  is  the  mean  of  the 
duplicate  samples  for  each  position. 
Difference  between  duplicates  was 
comparable  to  the  overall  difference 
observed,  suggesting  that  post-sample  optical 
geometry  is  not  a  significant  problem  for  the 
general  application  of  the  QFT. 
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Figure  2.  Analysis  of  the  relationship 
between  fi  and  ODf<A)  for  quadratic 
equation  curve  fits  including  different 
subsets  of  the  total  data  base  for  GF/F 
filters.  Exclusion  of  high  ODf<X)  values 
resulted  in  a  steeper  relationship,  while 
inclusion  of  all  data  up  to  ODf(X)  >  0.6 
caused  the  relationship  to  be  quite  flat.  An 
analysis  of  the  mean  determined  for 
discrete  ranges  is  included  as  the  discrete 
points  in  the  figure  where  fi  >  ODf<A)  / 
ODj(X).  The  ranges  for  discrete  analyses 
included  OD^fX)  from  O.O-O.OS;  O.OS-O.I; 
0.1-0.15;  0.15-0.2,  0.2-0.4;  and  0.4-0.6. 

The  points  are  plotted  in  the  middle  of  the 
range.  The  curve  fit  including  data  in  the 
range  0.0-0.4  best  fit  the  discrete  analysis. 


QFT  was  the  precision  of  the  technique  for  different  sized  particles,  or  for  particles  of  different 
refractive  index.  Although  it  has  been  hypothesized^*^  that  the  scattering  by  the  filter  itself  is  so 
intense  that  variations  in  particle  size  or  refractive  index  would  cause  insignificant  differences,  this  was 
not  previously  tested. 

A  detailed  study  of  the  QFT  for  GF/F  and  GF/C  filters  was  undertaken  with  a  variety  of 
phytoplankton  from  different  taxonomic  groups  and  of  different  cell  sizes.  The  phytoplankton  used  in 
the  analysis  included  Dunaliella  tertiolecta,  CMorella  sp.,  Nanochloroxis  sp.,  Thalassiosira  fluviatilis 
and  a  unicellular  cyanobacterium  Synechococcus  sp.  Among  the  phytoplankton,  the  diatom  is  expected 
to  have  the  highest  refractive  index  due  to  the  silicate  frustule.  Polystyrene  beads  (2.85  ^m  diameter) 
were  included  in  the  analysis  to  verify  the  validity  of  the  approach  for  high  refractive  index  particles. 

A  summary  of  the  effects  of  curve  fitting  data  covering  different  ranges  of  OD  is  illustrated  in 
Figure  2.  Values  of  ^  as  a  function  of  OD{<X)  are  presented  for  curve  fits  using  all  GF/F  datt  in  the 
ODf<X)  ranges  of  0-0.05;  0-0.4;  and  0-0.6.  Curve  fits  using  data  including  higher  values  of  OD^X) 
lead  to  a  flatter  response  for  p  due  to  the  lower  0  at  high  values  of  OD{<X).  Thus  the  low  OD  region 
was  compromised  for  a  better  fit  at  the  high  OD  end.  Conversely,  using  the  limited  lower  OD  range 
leads  to  fits  which  are  not  adequate  for  the  data  at  high  OD.  A  separate  analysis  of  0  determined  from 
the  experimental  observations  rather  than  the  curve  fit  is  also  presented.  The  discrete  points  plotted  in 
Figure  2  correspond  to  the  mean  value  of  the  ratio  of  OD^X)  /  OD,(X)  over  discrete  ranges  of  ODf<X), 
plotted  at  the  center  of  the  range.  It  is  evident  that  0  derived  from  the  curve  fit  data  using  the  range 
of  0-0.4  best  agrees  with  the  discrete  ratio  approach  up  to  a  maximum  OD  of  0.4.  This  is  the  basis  for 
restricting  the  final  curve  fitting  to  this  range.  From  a  practical  penpective,  one  generally  does  not 
work  with  samples  with  ODf(X)  >  0.4.  This  is  particularly  true  in  field  situations  where  particle 
concentrations  are  low  making  it  difficult  to  obtain  and  filter  sufficient  water  to  exceed  OD^X)  of  0.4. 
Thus  restricting  the  curve  fit  as  done  here  is  expected  to  yield  coefficients  of  the  most  general  utility. 
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GF/F 

GF/C 

n  observe 

lions 

6981 

5968 

a 

95%  lower 

0.380 

0.449 

estimate 

0.392 

0.459 

0.404 

0.469 

b 

95%  lower 

0.610 

0.602 

estimate 

0.655 

0.637 

0.700 

0.672 

Table  I.  Coefficients  for  the 
quadratic  equation  fit  (Equation 
9)  of  the  deuiled  studies  of  the 
relationship  between  ODj(A)  and 
ODf(A)  for  Whatman  GF/C  and 
GF/F  filter  types.  The  best 
estimate  and  the  9S%  confidence 
interval  for  the  coefficients  a 
and  b  are  provided. 


Data  from  this  analysis,  together  with  the  quadratic  equation  best  fit  are  presented  in  Figure  3a 
(GF/F  filter  type)  and  3b  (GF/C  filter  type).  Only  data  in  the  range  of  ODf(X)>-  0-0.4  were  used  in 
the  final  curve  fitting.  Coefficients  for  the  quadratic  equation  fit,  a  and  b,  and  the  95%  confidence 
limits  for  the  best  fit  are  provided  in  Table  1.  The  best  fit  results  for  GF/F  and  GF/C  are  presented 
in  Figure  4.  Figure  4a  illustrates  the  relationship  between  ODj(X)  and  ODf<X)  for  the  two  filter  types, 
and  Figure  4b  illustrates  the  relationship  between  and  ODf(A)  calculated  from  the  best  fit.  As 
reported  previously^,  0  for  GF/F  is  larger  than  0  for  GF/C,  and  both  show  a  non-linear  relationship 
with  ODf(X)  with  higher  values  at  low  ODf(X).  The  upper  limit  of  0  using  the  curve  fit  method  agrees 
with  the  upper  limit  of  ^  -  2.4  specified  previously ^  as  a  requirement  for  optimization  of  the  method 
so  that  over-correction  did  not  occur  at  low  values  of  ODf(X). 
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Figure  4.  (a)  A  comparison  of  the  best  fit  curves  for  the  relationship  between  ODj(A)  and  ODf(A)  for 
GF/C  and  GF/F  filters  (same  fit  as  3(a)  and  3(b).  (b)  A  comparison  of  the  relationship  between  /9, 
calculated  from  the  best  fit  in  (a),  and  ODj<X)  for  GF/C  and  GF/F  filter  types. 


2.3  Comparison  with  previous  GF/F  algorithm 

A  significant  body  of  literature  exists  using  the  correction  procedures  outlined  in  earlier  work. 
The  previous  procedure  for  the  commonly  used  GF/F  filters  was  indirect^’*®:  ODj<A)  for  the  GF/F 
was  converted  to  the  equivalent  ODf(A)  for  a  GF/C,  and  an  optimized  p  algorithm  for  GF/C  was 
applied  to  estimate  a  (A)  according  to  Equation  8.  The  present  procedure  provides  an  algorithm  for 
direct  estimate  of  ODj(A)  from  measurement  of  ODf<A)  for  various  filter  types. 

Comparisons  of  a  (A)  determined  from  measurement  of  ODj(A)  using  the  integrating  sphere,  and 
estimates  using  the  new  direct  algorithm  for  GF/F  or  the  previous  indirect  procedure  are  provided  for 
a  culture  and  a  sea  water  sample  in  Figure  Sa  and  Sb,  respectively.  The  sea  water  sample  was 
concentrated  by  filtration  and  resuspended  for  determination  of  ODj(A).  Small  volumes  of  the 
concentrated  suspension  (1-10  ml)  were  filtered  onto  GF/F  filters  for  ODf<A)  measurements. 

Estimates  of  ap(A)  using  the  previous  algorithm  are  as  low  as  75%  of  the  estimates  using  the 
new  algorithm.  Analysis  of  the  source  of  the  variability  indicates  that  the  relationship  previously  used 
to  estimate  ODf(A)  for  GF/C  from  ODf(A)  GF/F  is  responsible  for  most  of  difference.  Previously^’ 
linear  regressions  of  ODf<A)  for  GF/F  as  the  dependent  variable  and  ODf(A)  for  GF/C  as  the 
independent  variable  resulted  in  a  slope  of  1.3.  In  the  present  study,  the  slope  was  1.14.  The  previous 
algorithm  for  GF/C  provides  estimates  which  are  approximately  90%  of  the  estimate  of  the  new 
algorithm  for  GF/C,  which  is  within  the  precision  of  the  QFT.  Thus  it  appears  that  variation  between 
the  algorithms  is  due  to  the  conversion  from  GF/F  to  GF/C  using  the  previous  indirect  method.  It 
appears  that  different  filter  lots  cause  variability  in  the  algorithms  of  ±  15%. 
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Figure  S.  Comparison  of  estimates  of  a.CA)  using  the  integrating  sphere,  the  recommended  algorithm  ± 
the  95%  confidence  interval  for  the  GF/F  filter  type,  and  the  previous,  indirect  method  for  GF/F.  (a) 
Results  for  a  culture  of  Thalassiosira  fluviatilis.  (b)  Results  for  concentrated  sea  water  collected  from 
Scripps  Institution  of  Oceanography  pier.  Results  of  the  previous  algorithm  (lowest  curve  in  each  plot) 
underestimates  ap(A)  by  15-25%  compared  to  either  the  current  algorithm  or  the  integrating  sphere 
observations  on  suspensions. 


This  magnitude  of  variability  between  algorithms  is  greater  than  the  precision  determined  in  the 
analysis  of  variability  of  replicates  of  ODf(A)  measured  with  different  post-sample  acceptance  angles 
(Figure  1).  Volume  measurement  errors  are  expected  to  be  approximately  5%  and  errors  in  distribution 
of  particles  on  the  filter  are  probably  of  similar  magnitude.  Given  these  uncertainties,  the  QFT 
appears  to  have  an  overall  precision  of  ±  10%.  Error  in  curve  fitting  relationships  between  OD^(X)  and 
ODf(A)  on  specific  lots  of  filters  results  in  an  additional  5%  error.  Thus,  the  overall  accuracy  of  the 
procedure  is  of  the  order  ±  1 5  provided  the  curve  fit  is  for  the  filter  lot  used.  However,  an  additional 
error  of  comparable  magnitude  may  be  caused  by  variations  between  filter  lots,  and  if  this  is  not 
considered,  the  aggregate  error  may  be  as  great  as  20-30  %. 

2.4  Variations  among  fiiters  types 

It  is  not  possible  to  use  a  single  algorithm  for  the  QFT  for  samples  concentrated  on  different 
filter  types.  Differences  among  filter  types  have  been  reported^*^^  and  were  attributed  to  variations  in 
the  intensity  of  scattering  due  to  filter  characteristics  (thickness,  density  of  weave,  diameter  of  fibers, 
etc.).  A  direct  comparison  between  filter  types,  using  cultures  of  D.  tertiolecta,  CMorella  sp.  and  T. 
fluviatilis  was  carried  out.  The  data  used  for  this  analysis  are  a  subset  of  the  total  data  set  for  which 
the  same  volume  of  the  cultures  were  determined  for  all  of  the  four  filter  types.  For  each  species, 
ODj(A)  was  determined  on  a  dilution  of  the  culture  so  that  the  maximum  ODj(A)  in  a  1  cm  cuvette  was 
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between  0.05  and  0.1.  Several  volumes  covering  a  10-fold  range  were  filtered  onto  the  different  filter 
types.  The  original  ODj(A)  determined  in  a  1  cm  cuvette  was  scaled  by  /g  for  the  various  volumes 
filtered.  A  quadratic  equation  fit  to  the  data  was  made  for  each  filter  type  and  a  summary  of  the 
results  are  presented  in  Figure  6.  In  Figure  6a,  the  best  fit  curve  of  ODj(A)  as  a  function  of  OD{<A) 
for  each  filter  type  is  shown.  Figure  6b  demonstrates  that  ^  is  a  non-linear  function  of  OD^^A)  and  is 
different  for  the  various  filter  types.  The  best  fit  coefficients  for  this  analysis  are  presented  in  Table 
2.  The  best  fit  coefficients  for  GF/C  and  GF/F  are  slightly  different  than  those  for  the  larger  data  set 
presented  in  Table  I.  However,  in  the  range  of  ODf(A)  from  0.05  to  0.4,  the  two  fits  predict  ODj(A) 
within  10%  which  is  considered  to  be  the  overall  precision  of  the  QFT. 

Highest  values  of  P  were  observed  for  the  934AH  filter  type.  This  is  a  surprising  result  and 
does  not  agree  with  previous  observations^**^.  This  may  be  attributable  to  a  lower  tolerance  for  its 
manufacturing  specification  compared  to  the  other  filter  typbs.  The  allowable  manufacturing  tolerance 
for  range  in  thickness,  as  a  fraction  of  the  mean  thickness,  is  0.46,  0.26,  and  0.27  for  the  934AH, 
GF/C,  and  GF/F  filters  types,  respectively  (E.  Heilweil,  Whatman,  Inc.  technical  service,  personal 
communication).  The  magnitude  of  0  is  expected  to  increase  with  increased  thickness  of  the  filter. 
Previously,  the  Millcpore  HA  filter  was  determined  to  have  the  highest  0.  Now,  the  HA  type  is 
between  the  GF/F  and  the  GF/C.  As  previously  discussed^**^^  Ihg  filter  type  is  difficult  to  use 
because  it  is  not  as  hygroscopic  as  the  glass  filter  types.  Thus  it  is  difficult  to  keep  saturated  with 
water  unless  a  portion  is  actually  immersed  in  fluid.  This  was  not  possible  using  the  mounting 
procedure  previously  used,  but  the  present  method  did  maintain  saturation.  The  present  algorithm  for 
HA  is  therefore  considered  more  accurate,  provided  the  sample  and  reference  filters  are  kept  saturated 
during  the  analysis.  The  relationship  for  the  GF/C  filter  type  is  similar  to  previous  results. 


Figure  6.  Comparison  of  algorithm  results  for  Whatman  glass  fiber  filter  types  GF/C,  GF/F  and 
934AH  and  Millepore  cellulose  acetate  type  HA.  (a)  Curves  for  the  best  fit  of  the  quadratic  equation 
for  the  relationship  between  ODj(A)  and  ODf(A)  for  experiments  conducted  with  Dmaliella  tertiolecta, 
CMorella  sp.  and  Thalasiossira  fluviatilis.  Quadratic  equation  best  fit  coefficients  are  from  Table  2. 
(b)  The  relationship  between  0  and  ODj<A)  for  the  best  fit  relationships  presented  in  Figure  6a. 


SPIB  Vol  1 302  Ocean  Optics  X  (1 990)  /  14S 


GF/F 

GF/C 

934AH 

HA 

n  observa 

lions 

2800 

2800 

2800 

2800 

a 

95%  lower 

0.407 

0.491 

0.286 

0.430 

estimate 

0.415 

0.500 

0.298 

0.445 

0.422 

0.509 

0.309 

0.460 

b 

95%  lower 

0.665 

0.516 

0.663 

0.560 

estimate 

0.690 

0.547 

0.706 

0.611 

0.714 

0.577 

0.750 

0.663 

Table  2.  Coefficients  for  the  quadratic  equation  fit  (Equation  9)  of  the  limited  studies  of  the 
relationship  between  ODj(X)  and  ODf(Ji)  for  Whatman  GF/C,  GF/F,  934AH  and  Millepore  HA 
filter  types.  The  best  estimate  and  the  95%  confidence  interval  for  the  coefficients  a  and  b  are 
provided. 


2.5  Recommended  procedure 

Unlike  the  term  $  in  Equation  8,  the  coefficients  a  and  b  in  Equation  9  have  no  physical 
significance.  Nevertheless,  predicting  ODj(X)  from  ODf(A)  using  the  quadratic  equation  approach  is 
more  direct.  Furthermore,  estimates  of  are  subject  to  compounded  errors  by  taking  the  ratio  of 
experimentally  determined  values  of  ODj(A)  and  ODj<A)  in  the  algorithm  development.  Thus,  when 
using  Equation  8  as  previously  recommended  ••2, 16  estimates  of  ODj(A)  will  be  subject  to  uncertainty  in 
both  the  experimentally  determined  ODf(A),  and  the  compounded  error  in  estimates  of  Errors  in  $ 
are  of  greatest  concern  at  low  values  of  ODf(A)  where  instrumentation  noise  and  baseline  correction 
errors  become  a  more  significant  part  of  the  signal.  To  minimize  error  in  the  curve  fitting,  it  is  better 
to  fit  the  relationship  between  OD,(A)  and  OD^A)  rather  than  the  relationship  between  0 
(OD^A)/ODj(A))  and  ODj(A).  Once  a  fit  is  obtained,  it  can  be  solved  for  0,  if  desired  (Figures  2,  4b, 
6b)  .  It  is  therefore  recommended  that  estimates  of  ap(A)  be  based  on  a  direct  estimation  of  OD,(A) 
from  experimentally  determined  values  of  ODf(A)  using  the  quadratic  equation  fit  procedures  outlined 
above.  From  ODj(A),  ap(A)  can  be  calculated: 

ap(A)  -  2.3  0D5(A)  /  /g  (10) 

Some  practical  considerations  should  also  be  mentioned.  Unsatisfactory  filtration  results,  due  to 
clogging,  have  been  obtained  using  filtration  devices  with  a  fritted  glass  support  for  the  filter.  Better 
results  are  obtained  with  filtration  devices  which  have  either  stainless  steel  or  plastic  screens  supporting 
the  filter.  Once  filtered,  samples  should  be  kept  saturated  and  refrigerated  in  the  dark  until  analyzed. 

Previous  experiments^^  indicated  no  significant  loss  of  absorption  for  filtered  samples  of  a 
culture  kept  at  4  *C  for  up  to  24  hours,  but  changes  might  be  species  specific  depending  on  pigment 
degradation  enzymes  and  other  factors.  Significant  loss  of  absorption  (up  to  25%,  particularly  in  the 
red  peak)  and  spectral  shifts  in  the  blue  have  been  observed  for  filtered  samples  kept  frozen  for  6 
months  at  -20  *C  in  air  (B.  G.  Mitchell,  unpublished  data).  Until  a  satisfactory  procedure  can  be 
demonstrated  for  long-term  storage  of  samples,  it  is  strongly  recommended  that  analyses  be  done  as 
soon  as  possible  after  sampling. 
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Treatment  of  baseline  offsets  is  also  an  important  consideration.  In  general,  the  baseline  is 
spectrally  flat,  although  it  can  vary  by  ODf(A)  -  ±  0.03.  This  is  of  no  consequence  provided  a 
wavelength  with  no  absorption  is  available  as  a  null  reference.  For  cultures  such  as  those  studied  here, 
the  assumption  of  no  absorption  at  750  nm  is  reasonable.  However,  this  assumption  poses  problems  in 
field  work,  particularly  if  a  significant  fraction  of  the  particulates  are  of  detrital  origin.  Nevertheless, 
considering  baseline  shifts  of  ODf(A)  ±  0.03,  the  magnitude  of  error  if  a  750  nm  null  correction  is  not 
used  would  be  greater  than  the  errors  associated  with  implementation  of  the  usually  reasonable 
assumption  that  absorption  at  750  nm  is  zero. 

3.  CONCLUSIONS 

A  detailed  analysis  was  conducted  regarding  the  generality,  accuracy  and  precision  of  the 
quantitative  filter  technique  (QFT)  for  determining  the  absorption  coefficient  of  aquatic  particulates, 
ap(A).  The  generality  of  the  procedure  was  verified  through  the  observation  that  changing  the  angle  of 
acceptance  to  the  detector  of  the  post-sample  diffuse  light  did  not  significantly  alter  the  magnitude  of 
ODf(A),  provided  the  baseline  was  properly  determined.  The  QFT  is  capable  of  estimating  the 
magnitude  of  ap(A)  with  a  precision  of  better  than  10%  and  an  overall  accuracy  within  t  IS  %.  This 
accuracy  is  attained  if  coefficients  for  a  specified  lot  of  filters  are  determined.  The  non-linearity 
between  p  and  ODf(A)  was  verified  as  was  the  difference  in  $  for  different  filter  types. 

A  new,  more  direct,  algorithm  for  the  GF/F  filter  type  results  in  15-25%  higher  estimates  of 
ap(A)  compared  to  the  previous  indirect  method  of  converting  ODf(A)  for  GF/F  to  ODf(A)  for  GF/C 
and  then  correcting  with  p  for  the  GF/C.  Although  the  GF/C  and  GF/F  filter  types  appear  to  be 
more  consistent  with  previous  observations,  a  difference  in  the  slope  of  linear  fits  between  ODf<A)  for 
GF/F  and  GF/C,  compared  to  previous  observations^,  is  hypothesized  to  be  caused  by  variations  in 
filter  lots.  The  higher  slope  previously  observed  results  in  an  overcorrection  when  applied  to  analyses 
with  the  current  stock  of  filters.  Results  in  this  study  for  the  934AH  filter  type  were  significantly 
different  than  previous  results^M^.  Apparently,  filters  from  different  lots  have  different  thickness 
which  may  cause  different  relationships  for  p.  If  these  differences  are  not  considered,  an  additional 
error  of  10-15  %  beyond  the  accuracy  of  the  procedure  may  result.  For  the  highest  possible  accuracy 
of  the  QFT,  it  may  be  necessary  to  determine  the  coefficients  for  each  lot  of  filters  used. 
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